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A Particle Multi-Target Tracker for Superpositional 
Measurements using Labeled Random Finite Sets 

Francesco Papi and Du Yong Kim 


Abstract —In this paper we present a general solution for multi¬ 
target tracking with superpositional measurements. Measurements 
that are functions of the sum of the contributions of the 
targets present in the surveillance area are called superpositional 
measurements. We base our modelling on Labeled Random Finite 
Set (RFS) in order to jointly estimate the number of targets and 
their trajectories. This modelling leads to a labeled version of 
Mahler’s multi-target Bayes filter. However, a straightforward 
implementation of this tracker using Sequential Monte Carlo 
(SMC) methods is not feasible due to the difficulties of sampling 
in high dimensional spaces. We propose an efficient multi-target 
sampling strategy based on Superpositional Approximate CPHD 
(SA-CPHD) filter and the recently introduced Labeled Multi- 
Bernoulli (LMB) and Vo-Yo densities. The applicability of the 
proposed approach is verified through simulation in a challenging 
radar application with closely spaced targets and low signal-to- 
noise ratio. 

Index Terms —Labeled RFS, Superpositional Measurements, 
CPHD Filtering, Proposal Distribution 


I. INTRODUCTION 

Superpositional sensors are an important class of pre¬ 
detection sensor models which arise in a wide range of joint 
detection and estimation problems. For example, in problems 
such as direction-of-arrival estimation for linear antenna arrays 
(1 j, multi-user detection for wireless communication networks 
(2j, acoustic amplitude sensors (3), radio frequency (RF) 
tomography 0, target tracking with unresolved or merged 
measurements (5) , multi-target Track-Before-Detect with 
closely spaced targets 0- 0- the sensor output is a func¬ 
tion of the sum of contributions from individual sources. In 
classical estimation theory, a frequency domain model for the 
superpositional sensor is generally used to design algorithms 
for source separation and parameter estimation. Conversely 
in dynamic state estimation, a detection based model is typ¬ 
ically employed to transform the collected data into a set of 
point measurements, in order to facilitate the development 
of computationally efficient estimation algorithms. This is 
specifically the case in multi-target tracking |9)-]TT|, which is 
an important problem in estimation theory involving the joint 
estimation of an unknown and time varying number of targets 
and their trajectories. 
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Many real life applications in radar, sonar 0 -PI, com¬ 
puter vision [f3||-p3|, robotics G9-G3. automotive safety 
0, 0, cell biology (22j-(26|, etc., can be described 
as multi-target tracking problems. Most of the multi-target 
tracking algorithms existing in the literature are designed for 
data that have been preprocessed into point measurements 
or detections 0-0, 0. These algorithms are based on 
the “detection sensor” model which assumes that each target 
generates at most one detection, and that each measurement 
belongs to at most one target (TO). The performed preprocess¬ 
ing of raw measurements into a finite sets of points is efficient 
in terms of memory and computational requirements, and is 
usually effective for a wide range of applications. However, 
the compression might lead to significant information loss in 
the presence of low signal-to-noise ratio (SNR) and/or closely 
spaced targets. The standard “detection sensor” approach may 
not be adequate in this case, and making use of all information 
contained in the pre-detection measurements becomes neces¬ 
sary. In turn, this requires more advanced sensor models and 
new algorithms. 

In a superpositional sensor model, the measurement at 
each time step is a superposition of measurements generated 
by each of the targets present |28|. In (28) Mahler derived 
a superpositional Cardinalized Probability Hypthesis Density 
(CPHD) filter as a tractable approximation to the Bayes multi¬ 
target filter for superpositional sensor. The approach was im¬ 
plemented in [4) using SMC methods, and successfully applied 
to a passive acoustics application as well as RF tomography. 
The technique was also extended to multi-Bernoulli and a 
combination of multi-Bernouli and CPHD |[29j, (30) . These 
filters, however, are not multi-target trackers because they rest 
on the premise that targets are indistinguishable. Moreover, 
they require at least two levels of approximations: analytic 
approximations of the Bayes multi-target filter and particles 
approximation of the obtained recursion. 

Inspired by 0, 0, this paper proposes a multi-target 
tracker for superpositional sensors which estimates target 
tracks and requires only one level of approximation. Our 
formulation is based on the same random finite set (RFS) 
framework that the superpositional CPHD filters 0, (28) were 
derived from. However, we used a special class of RFS models, 
called labelled RFS 0- which enables the estimation of 
target tracks as well as direct particle approximation of the 
(labeled) Bayes multi-target filter. To mitigate the depletion 
problem arising from sampling in high dimensional space we 
propose an efficient multi-target sampling strategy using the 
superpositional CPHD filter 0. In particular, we will show 
how the recently introduced Labeled Multi-Bernoulli (32) and 
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Vo-Vo Q 1341 densities can be constructed from the superposi- 
tional approximate CPHD (SA-CPHD) filter. These densities 
are then used to design effective proposal distributions for 
the RTS multi-target particle filter. While both the CPHD and 
labeled RFS solutions require particle approximation, the latter 
has the advantage that it does not require particle clustering 
for the multi-target state estimation. The applicability of the 
proposed approach is verified through simulation analyses in 
a challenging closely-spaced multi-target scenario using radar 
power measurements with low signal-to-noise ratio (SNR) 


The paper is organised as follows: in Section [II] we re¬ 
call some definitions for Labeled RFSs and superpositional 
sensors. In Section III we discuss the multi-target particle 
filter, the labeled multi-target transition density, and the su¬ 
perpositional approximate CPHD. Two multi-target particle 
trackers using the Labeled Multi-Bernoulli (LMB) and the 
Vo-Vo densities for the proposal distribution are presented 
in Section IV Numerical results for a radar application are 
presented in Section [VJ while conclusions and future research 
directions are discussed in Section PvTl 


II. Background 


This section briefly presents background material on su¬ 
perpositional sensor model and the RFS framework which 
we adopt for the formulation of a multi-target tracking filter. 
Subsection |II-A| provides a summary of basic concepts in RFS. 
We present a concise description of the superpositional sensor 
model in subsection II-B and report a summary of key ideas 
on labled RFS needed for the derivation in subsection 1II-CI 


A. Multi-target Estimation 

Suppose that at time k, there are N(k) target states 
Xk,i, ■ ■ ■, Xk,N(k)> eac h taking values in a state space X. 
In the random finite set (RFS) framework, the multi-target 
state at time k is represented by the finite set X k = 
{xk,i, ■ ■ ■, Xk,N(k)}' and the multi-target state space is the 
space of all finite subsets of X, denoted as T(X). An RFS 
is simply a random variable that take values the space T(X~) 
that does not inherit the usual Euclidean notion of integration 
and density. Mahler’s Finite Set Statistics (FISST) provides 
powerful yet practical mathematical tools for dealing with 
RFSs (Tg, |36) based on a notion of integration/density that 
is consistent with point process theory |37j. 

Similar to the standard state space model, the multi-target 
system model can be specified, for each time step k, via 
the multi-target transition density fk\k-i and the multi-target 
likelihood function g k , using the FISST notion of integra¬ 
tion/density. The multi-target posterior density (or simply 
multi-target posterior) contains all information about the multi¬ 
target states given the measurement history. The multi-target 
posterior recursion is direct generalisation of the standard 

'The Vo-Vo density was originally called the Generalized Labeled Multi- 
Bernoulli density. However, for compactness we follow Mahler’s latest book 
[33[ and call this the Vo-Vo density 


posterior recursion {38), i.e. 

Tto-.k{X 0 . k \z X:k ) oc (1) 

gk(Zk\Xk)fk\k-l(Xk\Xk-l)TtO:k-l(X 0: k-l\Zi : k-l) 

for k > 1, where X 0:k = (X 0 , ..., X k ), and z\, k = (zi ,..., z k ) 
is the measurement history with z k denoting the measure¬ 
ment vector at time k. Target trajectories or tracks can be 
accommodated in the RFS formulation by incorporating a 
label in the target’s state vector (10), (3T|, (39). The multi¬ 
target posterior 0 then contains all information on the random 
finite set of tracks, given the measurement history. In [39|, 
the set of tracks are estimated by simulating from the multi¬ 
target posterior 0 using particle Markov Chain Monte Carlo 
(PMCMC) techniques (40) . 

Computing the multi-target posterior is prohibitively expen¬ 
sive for on-line applications. A more tractable alternative is 
the marginal n k at time k known as the multi-target filtering 
density. For notational compactness we omit the dependence 
on the measurement history. Marginalizing the multi-target 
posterior recursion 0 yields the multi-target Bayes filter (lO) , 

HD’ 

7T k \ k -i(X k ) = J f klk _ 1 (X k \X)7T k . 1 (X)6X, ( 2 ) 

/ y \ _ 9k(z k \X k )Tt k \ k -i{X k ) 

nk[ k) ~7 ;,i 

where it k \k-i is the multi-target prediction density to time 
k, and the integral is a set integral defined for any function 
f:F{X) -» «. by 

« OO 

/ f(X)SX = 

2 = 0 

In (31), (34) an analytic solution to the multi-target Bayes filter 
0 , 0, known as the Vo-Vo filter (33), was derived using 
labeled RFSs. Note that the majority of work in multi-target 
tracking is based on filtering, and often the term "multi-target 
posterior" is used in place of "multi-target filtering density". 


Tj J f({x i,..., Xi})d(x i,..., xi). 


(4) 


B. Superpositional Sensor 

In a superpositional sensor model, the measurement z is 
a non-linear function of the sum of the contributions of 
individual targets and noise, i.e. 

z = $ H 3'( x )’ e j (5) 

Wx / 


where j(x) represents the contribution of the single-target 
state x to the sensor measurement (7 is a non-linear mapping 
in general), e is the measurement noise, and $ is a nonlinear 
mapping. For example, a superpositional sensor model com¬ 
monly used in radar is 


z = 


e je ^2 A(x)g(x) + C 

x&X 


2 


( 6 ) 


where g(x) is the point-spread function of target x, A(x) is the 
(known) amplitude, 9 is the phase noise, uniformly distributed 
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on [0, 2w], and £ is circularly complex symmetric Gaussian 
noise. It is clear that this model takes on the form 0 by 
defining e = (0,£). In general, the multi-target likelihood 
function for the superpositional sensor model is the probability 
density of the measurement z given YhxeX the sum °f 
the contributions of individual targets, i.e. 

gk(z\X) = h k (*; ) ’ (7) 

V xex ) 

The SA-CPHD filter filter presented in ||4j is an approximation 
to the multi-target Bayes filter for a superpositional measure¬ 
ment model of the form 

z= 'Yh T'fc) + e ’ (8) 

xex 

where e is distributed according to Mr, a zero mean Gaussian 
with covariance R. Hence, the likelihood function for the 
superpositional measurement is 

9 k{z\X) = Mr f z- ^2 j(x) I . (9) 

V xex ) 

Similar to the CPHD filter (for the standard sensor model) 
|4T) , the SA-CPHD filter filter [4j is an analytic approximation 
of the Bayes multi-target filter 0,0 based on independently 
and identically distributed (iid) cluster RFS. A brief review of 
the SA-CPHD filter is given in subsection |III-C| Both filters 
recursively propagate the cardinality distribution and the PHD 
of the posterior multi-target RFS. The CPHD filter can be im¬ 
plemented with Gaussian mixtures or particles {42), while only 
the particle implementation is available for the SA-CPHD filter 
|4). Particle implementations of PHD/CPHD filter in general 
require clustering to extract multi-target estimates, which can 
introduce additional errors under challenging scenarios. 


C. Labeled RFS 


To perform tracking in the RFS framework we use the 
labeled RFS model which incorporates a unique label in the 
target’s state vector to identify its trajectory ©■ In this model, 
the single-target state space X is a Cartesian product XxL, 
where X is the feature/kinematic space and L is the (discrete) 
label space. A finite subset set X of XxL has distinct labels 
if and only if X and its labels {£: [x, £) £ X} have the same 
cardinality. An RFS on XxL with distinct labels is called a 
labeled RFS 0D, (34). 

For the rest of the paper, we use the standard inner product 
notation (f,g) = J f(x)g(x)dx. We denote a generalization 
of the Kroneker delta and the inclusion function that take 
arbitrary arguments such as sets, vectors, by 



1, if X = Y 
0 , otherwise 



1, if X C Y 
0 , otherwise 


We also write ly(x) in place of ly({x}) when X = {a;}. 
Single-target states are represented by lowercase letters, e.g. x, 
x while multi-target states are represented by uppercase letters, 
e.g. A', X, symbols for labeled states and their distributions 
are bolded to distinguish them from unlabeled ones, e.g. x. 


X, 7 r, etc, spaces are represented by blackboard bold e.g. X, 
Z, L, etc. 

An important class of labeled RFS distribution is the 
generalized labeled multi-Bernoulli distribution (31 i), known 
as the Vo-Vo distribution (331, which is the basis of an 


analytic solution to the Bayes multi-target filter 1341. Un¬ 
der the standard multi-target measurement model, the Vo- 
Vo distribution is a conjugate prior that is also closed under 
the Chapman-Kolmogorov equation. If we start with a Vo- 
Vo initial prior, then the multi-target posterior at any time is a 
also a Vo-Vo distribution. Let C : XxL — > L be the projection 
C{{x,£)) = £, let A(X) = <5| X | (|jC(X)|) denote the distinct 

label indicator , and h x = h(x), denote the multi-object 

xex 

exponential, where h is a real-valued function, with \v = 1 
by convention. A Vo-Vo density is a labeled RFS density on 

J'(XxL) 


7r(X) = A(X)^ w ( c )(£(X)) 
cec 


,(c) 


n x 


GO) 


where C is a discrete index set, w^(L) and p^ satisfy: 

EE w(c) ( L ) = 1 , ( 11 ) 

-LCLceC 

/jXW)* = 1. (12) 

The Vo-Vo density (TO) can be interpreted as a mixture of 
multi-object exponentials. Each term in © consists of a 
weight u/ c )(£(X)) that depends only on the labels of X, and 
a multi-object exponential [p <c> ] that depends on the entire 
X. The Labeled Multi-Bernoulli (LMB) family is a special 
case of the Vo-Vo density with one term of the form: 


7T(X) 

= A(X)tu(£(X))p x 

(13) 

p(x,£) 

II 

1 ? 

(14) 

uj{L) 

= n on tS? 

reM teL 

(15) 


where {(r^,pW)}ee m, M C L, is a given set of parameters 
with representing the existence probability of track £, and 
p^> the probability density of the kinematic state of track £ 
given its existence (3T). Note that for an LMB the index space 
C has only one element, in which case the (c) superscript is 
not needed. The LMB family is the basis of the LMB filter, 
a principled and efficient approximation of the Bayes multi¬ 
target tracking filter, which is highly parallelizable and capable 
of tracking large number of targets (43) . 


III. Bayesian multi-target tracking for 

SUPERPOSITIONAL SENSOR 

In this section we describe the classical particle Bayes multi¬ 
target filter which has very high computational com¬ 

plexity in general. Fortunately, using labeled targets greatly 
simplifies the multi-target transition density and drastically re¬ 
duces the computational complexity. Subsection |III-A| presents 
a summary of the classical multi-target particle filter and 


Subsection III-B details the labeled multi-target transition 
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density that reduces the computational complexity. Subsection 
|III-C| reviews the equations of the superpositional CPHD filter 
that is used to following section to construct LMB/Vo-Vo 
efficient proposal distribution for the multi-target particle filter. 

Following 0,0, to ensure distinct labels we assign each 
target an ordered pair of integers £ = (k, i), where k is the time 
of birth and i is a unique index to distinguish targets born at the 
same time. The label space for targets born at time k is denoted 
as Lfc, and a target born at time k, has state x € XxLfc. The 
label space for targets at time k (including those born prior 
to k), denoted as Lo : fc, is constructed recursively by Lo : fc = 
Lo;fc_i U Lfc (note that Lo : fc_i and L/, are disjoint). A multi¬ 
target state X at time k, is a finite subset of X = XxLo : fc. 
For completeness, the Bayes multi-target tracking filter, i.e. 
the multi-target Bayes recursion ([2]), Q for labeled RFS, is 
provided below 

*fc|fc-t(Xfc) = J ^| fc _ 1 (Xfe|X)7Tfe_ 1 (X)(5X, 

/.y. \ Qk^Zkl'X-k)'^k\k— l(Xfc) 

k = f g k (z k \X)7r k , k _ 1 (X)SX- 

A. Particle Bayes multi-target filter 

The propagation of the multi-target posterior involves the 
evaluation of multiple set integrals and hence the compu¬ 
tational requirement is much more intensive than single¬ 
target filtering. Particle filtering techniques permit recursive 
propagation of the set of weighted particles that approximate 
the posterior. Central in Monte Carlo methods is the notion of 
approximating the integrals of interest using random samples. 
While the FISST density is not a density (in the Radon- 
Nikodym context), it can be converted into a probability 
density (with respect to a particular dominating measure) by 
cancelling out the unit of measurement 0. Monte Carlo 
approximations of the integrals of interest can then be con¬ 
structed using random samples. The single-target particle filter 
can thus be directly generalised to the multi-target case. In the 
multi-target context however, each particle is a finite set and 
the particles themselves can thus be of varying dimensions. 
Following (37) , suppose that at time k — 1, a set of weighted 
particles xi-i}^ representing the multi-target pos¬ 

terior 7rfc_ 1 |fe_ 1 is available, i.e. 

N p 

Tr^Xj^^^XjX^) (18) 

i=1 

Note that i5(-;X^ i l 1 ) is the Dirac-delta concentrated at xj^r 
(different from the Kronecker-delta Sx that takes values of 
either 1 or 0). The particle filter proceeds to approximate the 
multi-target posterior at time A: by a new set of weighted 
particles {ufi > , xj c , ' , }^j as follows 

Multi-target Particle Filter 


(16) 

(17) 


For time k > 1 


For i = 1.. .., N p sample xjj* ~ q(-|Xj*j_ ] , zy.) and set 
_ 9fe ( 2fc |x6 ) )f fe | fe _ 1 (x6 ) |X« 1 ) 




-(») 

w k 


q fe (X«| Xj^.Zfc) 


M 

k -1 


(19) 


• Normalize the weights: w) ’ — 




E 2 




Resampling Step 

• Resample to S et , X^}^ 

The importance sampling density qfc(-|Xfc_ 1 , zy.) is a multi¬ 
target density and X/, is a sample from an RFS. It is implicit 
in the above algorithm description that 

ffclfc-itXfcIXfc-ij 


sup 

X fc ,X fc _! 


Hfc (X/-|Xfc_i, ~fc) 


< 00 


( 20 ) 


so that the weights are well-defined. Convergence results for 
the multi-target particle filter are given in (37) . 

Notice that the entire posterior can be computed by mod¬ 
ifying the pseudo-code of the multi-target particle filter so 


that Xy' ; fc is used in place of Xjx 1 and X, 
in place of Xq.[,_ 1 . This would in principle solve the so 


(0 


r {i) 

0:fe —1 


is used 


called mixed labelling problem |44|. However, this is com¬ 


putationally demanding because it requires recomputing the 
whole history of each multi-target particle (39) . Alternatively, 
forward-backward smoothing can be used to approximate the 
entire posterior 0 . In this paper we focus on designing 
efficient proposal distributions for the multi-target particle 
filter approximating the filtering recursion. In future work we 
will consider the application of the proposed approach to the 
problem of estimating the full posterior. 

The main practical problem with the multi-target particle 
filter is the need to perform importance sampling in very 
high dimensional spaces if many targets are present. In (32) , 
|45|. |46J, the transition density is used as the proposal, 
i.e. qfc(-|X^l lf z fc ) = f fe | fc _i(jX^J. While this avoids the 
evaluation of the transition density, it suffers from particle 
depletion even for a small number of targets. This prob¬ 
lem is compounded with superpositional sensor due to less 
informative measurements arising from low SNR. A naive 
choice of importance density such as the transition density 
will typically lead to an algorithm whose efficiency decreases 
exponentially with the number of targets for a fixed number 
of particles (37). The problem with using a proposal other 
than the transition density, is that the weights are difficult 
to evaluate due to the combinatorial nature of the transition 
density for unlabled RFS. Fortunately, for labeled RFS the 
transition density simplifies to a form that is inexpensive to 
evaluate. 


B. Labeled multi-target transition density 

The multi-target transition model for labeled RFS is sum¬ 
marised as follows. Given a multi-target state X at time fc — 1, 
each state (x,£) £ X either continues to exist at the next 
time step with probability ps{x, £) and evolves to a new state 
(xfc, tk) with probability density f(xk\x, l)8i(£k), °r dies with 
probability 1 —ps(x, i). In addition, the set of new targets born 
at time k is distributed according to the LMB distribution 

bfc(Y) = A(Y)w B , fc (£(Y)) [p B , fc ] Y ( 21 ) 

where u)B,k a 'id p B k are given parameters of the multi-target 
birth density b^, defined on T(Xx L^). Note that bfc(Y) = 0 
if Y contains any element y with £(y) ^ L k- The birth model 















5 


( |2T| covers both labeled Poisson and labeled multi-Bernoulli 
[311. The multi-target state X&, at time k, is the superposition 
of surviving targets and new born targets. The model uses the 
standard assumption that targets evolve independently of each 
other and that births are independent of surviving targets. 

It was shown in [311 that the multi-target transition density 
is given by 

ffc|fc-t (Xfc|X) = 

s fc | fc _ 1 (Xfc (T (X x L 0 : fc-i)|X)bfc(Xfc (T (X x Lfc)) 

( 22 ) 


where 


Sfe|/b- 1 (W|X) 




A(W)A(X)1 £(x) (£(W)) [$ fc | fc _i(W; -)] X 

(23) 

/ ps(x,£) f(xk\x,£), if (xk,£) 6 W 
\ 1 -ps(x,£), if£#£(W) 

(24) 


Unlike the general multi-target transition density (see (B 
|), the special case for labeled RFS (22 1 does not contain 


any combinatorial sums. It is simply a product of terms 
corresponding to the surviving targets and new targets. Con¬ 
sequently, numerical complexity of the weight update in the 
multi-target particle filter drastically reduces. 


C. Superpositional Approximate CPHD filter 


In this section we recall the approximate CPHD for super- 
positional measurements of the following form: 


Zk 




2 

+ n fc 


(25) 


PHD, while the update step of the SA-CPHD filter is given 
by: 


^ , , A/- Sr+E ~ (z fe - nA fc | fe -i) 

Pk(n) =Pk\k-i(n)— r -t -—-1-v ( 26 ) 


Vk{x) =v k \ k -i(x) 


A/e r +S fc | fc _ 1 ^k\k— 1 frk\k— l) 

(z k -h{x)-p° k |fc-1 ) 
■^T lr +T lk \ k _ 1 (Zk Nk\k-lPk\k-l) 


(27) 


where X r = cr n Im is the noise covariance, N k \k-i is the 
predicted average number of targets and: 


Pk= J h(x)s k \ k -i{x)dx 
Sfc = J h(x)h(x) J s k \ k _ 1 (x)dx 
T, k =n ( e * - PkP-l'j 

Efc =N k \ k _ 1 tl k + — N k \ k _ 1 p k p k ) 


Pk 


Gj%-iW 

Nk\ k- 


k\k-l 


"A L k 


r '( 2 ) (i'i 

^/c|fc-l V- 1 -; - 

^ k ^ k + 


gfclLit 1 ) 

Nkik-i 


N k it_i 


(28) 

(29) 

(30) 

(31) 

(32) 


PkP k 

(33) 


where SjUfc_i(:r) is the normalized predicted intensity, and 
&l\ k - V Gfcifc_i(l) and '-'fcife-i ( 1 ) are th e variance, second 
factorial moment and third factorial moment of the predicted 
cardinality distribution Pk\k-i( n )- The equations of the su¬ 
perpositional approiximate CPHD filter can be implemented 
efficiently using SMC methods. In the following section we 
describe how the updated PHD and cardinality distribution 
from the SA-CPHD filter can be used to design efficient 
proposal distributions for multi-target tracking. 


where X& is the multi-target state at time k, n k ~ A/"(0, <r^) 
is zero-mean white Gaussian noise, and h(x) is a possibly 
nonlinear function of the single state vector x. Notice that the 
model in eq. ( |25j ) can be used to approximate the radar power 
measurement eq. © assuming a Gaussian noise in power. 
Obviously the model in eq. (251 is a strong approximation 
of eq. However, it allows using the update step of the 
SA-CPHD filter to evaluate measurement updated intensity 
function w(x) and cardinality distribution p(n) for the target 
set. In turn, the information in the updated u(x) and pin), 
along with the targets labels from the previous step and birth 
process, can be used to construct an approximate posterior 
density using the Vo-Vo and/or LMB distributions in eq. © 
and ©• Finally, the obtained approximate posterior is used 
as a proposal distribution for the multi-object particle filter. 

The vector measurement z k in eq. ( [25} usually represents an 
array of sensors for SA-CPHD filter, e.g. acoustic amplitude 
sensors, radio-frequency tomography, etc. For application of 
the SA-CPHD filter to tracking using radar power returns, 
the vector measurement z k contains the radar power returns 
from the set of cells being interrogated by the radar at time 


k. Hence, z k = [z^ ...z- 
being interrogated by the radar. Following ||4]|, standard CPHD 
formulas are used for the predicted cardinality distribution and 


(m)i 


where m is the number of cells 


IV. Efficient Proposal Distributions based on 
Superpositional Approximate CPHD filter 

In this section we detail the CPHD-based proposal distri¬ 
bution and the multi-target particle filter equations. In super- 
positional multi-target filtering, the multi-target posterior gen¬ 
erally cannot be written as a product of independent densities 
because the target states are statistically dependent through 
the measurement update. This means that an effective particle 
approximation of the posterior distribution is of great interest. 
Unfortunately, designing an effective multi-object proposal 
distribution q(Xfc|Xfc_i, z k ) is not a simple task when using 
superpositional sensors. In this section we exploit the SA- 
CPHD filter to construct a relatively inexpensive LMB based 
proposal as well as more accurate Vo-Vo based proposal. The 
basic idea is to obtain the updated PHD v k (x) and cardinality 
distribution p k (n) at time k from the SA-CPHD filter and 
construct a proposal distribution q(Xfc|Xfc_i, z k ) that exploits 
the approximate posterior information contained in both the 
cardinality distribution p k {-) and the state samples from v k (■). 

Assume a particle representation j X^-,, \ of the 

v J i—1 

posterior distribution 7 rfc_!(X) is available at time k— 1. Then, 
the cardinality distribution p(n) and the PHD v(x) of the 
unlabeled multi-target state at time k — 1 are given by J37): 
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Pk-i(n) oc 

X" (*) 

w k- 1 

(34) 


i: l x fe-ll= n 


Vk-l(x) = 

N p 

E E w k-i s ( x '’ x k-i,e) 

(35) 


1=1 


where x^}_ 1>e 

denotes the kinematic part of 

each 


(zk-v^k-i) e x[Zi, and is the Dirac-delta 

concentrated at x < A__ l v The superpositional CPHD is then used 
to obtain the update cardinality distribution Pk(n) and PHD 
Vk{x) using the measurement z k collected at time k. Notice 
that differently from standard unlabeled CPHD filtering, there 
is a natural labeling/clustering of particles due to the existing 
labels at time k — 1 and the chosen iid cluster process with 
implicit cluster labels for the birth model. In fact, let Vk{x) 
be the updated PHD at time k 

N p 

vk{x) = J2 w k )s ( x > 4‘i) < 36 ) 

i=1 fcc(x«) 

Then we can rewrite the (unlabeled) PHD as a sum over all 
labels Lo : fc of labeled PHD terms Vk,e{x), i.e. 

Vk{x) = ^2 v kA x ) (37) 

where 

N p 

Vk,e( x ) = E Se(£ ) w2 6 (x;xfy) (38) 

i=1 ^£(x«) 

is the contribution to the PHD of track £. Note that the 
above is not the PHD of a labeled RFS but the PHD mass 
from a specific label representing a survival or birth target. 
This means that at time k we can extract |Lo : fc| clusters of 
particles from the posterior PHD. Furthermore, a continuous 
approximation to each cluster can be obtained by evaluating 
sample mean and covariance for a Gaussian approximation 
to Vk,i(-)- Alternatively, it is possible to use kernel density 
estimation (KDE), however this will not be considered in this 
paper. For £ £ Lo : fc, let p k ,e and Qk,t denote the sample 
mean and covariance corresponding to the PHD cluster Vk,e (•)■ 
Hence, we approximate the PHD clusters as follows 


Vk,e ( x ) = Pk W ■A/’ (x-, p,k,i, Qk,e ) (39) 

where p ^ (£) is the PHD mass of the t th cluster. For the sake 
of explicitness, in our exposition we denote the PHD mass of 
survival targets as s (A and the PHD mass of newly born 


targets as Pfc B (f), 

Np 

Pfc, S W=E E iU£U:k-l (40) 

*=Ve,c(x«) 

N p 

Pfc,flW = E E ifieu (41) 

i=1 <'e£(x>‘>) 

In practice we constrain the survival and birth probabilities 

Pk,si ^ \PS,min\ PS,max ] and Pfc.sW e \p B ,min ; PB,max\- 

The constraint p. t min is imposed to avoid the complete loss of 
a track due to errors in the CPHD update while the constraint 
P-,max is required since the PHD in each track cluster can 
exceed 1. 

The obtained posterior cardinality and posterior target 
clusters can be used to construct a proposal distribution 
q(-|X*;_i, Zk). In the following subsections we detail two 
strategies for constructing the proposal q(-|Xfe_ 1 ,Zfe) as an 
LMB density of the form ( p~3] > and as a Vo-Vo density of the 
form ( flQ| . 


A. LMB Proposal Distributions 

In this subsection we describe how to construct a multi¬ 
target proposal distribution q (X/,.|X_-|, Zk) for the multi¬ 
target particle tracker by using an LMB density, i.e. 


q (Xfc|Xfc_i, Zk) — 

qs(X fc n (X x L 0: fc)|X fc _ 1 )q B (X fc — (X x Lo;fc)) (42) 


where q s (X k G (X x L 0:fe )|X fe _i) and q B (X fc - (X x L 0:fc )) 

are the LMB proposals for survival targets and birth tar¬ 
gets, respectively. Specifically, the survival and birth propos¬ 
als are constructed using the CPHD updated birth Pk,B W 
and survival p~£ s {£) probabilities and the Gaussian clusters 
p,k,e, Qk,t)- F° r the survival proposal we have. 


qs(Wfc|Xfc_i, Zk) = 

A(W fc )A(X fc _ 1 )l £(Xft _ l) (r(W fc )) [<MW fc ; •)] Xfc_1 (43) 


$s(Wk-,x,e) 


Pk sW N(x-,p.k,t,Qk,i), if (x,£) £ Wfc 
l 1 Pits W, ifitXWk) 

(44) 


while for the birth proposal we have. 


qs(Xfe) 


$s(X fc ) 


A (X fc )[$ B (.)] Xfc 
j Pk B (^) N{x', pk,l i Qk,t)i 

l i 1 p£ b (Q> 


(45) 

if (x,£) £ X k 
if £ i £(X k ) 

(46) 


In summary, the multi-target proposal distribution in eq. (19i 
is constructed using two LMB densities for the existing 
and newly appeared targets, respectively. A pseudo-code of 
the multi-target particle filter using the LMB proposal for 
sampling is given below. Notice that we used the following 
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definitions for grouping of labels in each particle 


r(i) 

= £(x« 1 )n^(x«) 

(47) 


= £(X« 1 )-£(X«) 

(48) 


= l,H£(x«) 

(49) 

r.(i) 

''NB 

= L fc -£(X«) 

(50) 


where for each particle ( i ), d'i' 1 is the set of survived labels, 
c!']] is the set of death labels, Cg is the set of labels for newly 
born targets, and C$ B is the set of labels that did not generate 
a new targets. 


a proposal distribution that is efficient in low SNR scenar¬ 
ios. Generally, for low SNR the Multi-Bernoulli cardinality 
distribution is not sufficiently informative, so that being able 
to estimate a more general cardinality distribution becomes 
fundamental. This reasoning is true also in classical multi¬ 
target tracking with detection measurements, e.g. the CPHD 
filter outperforms the Multi-Bernoulli filter in low SNR fTO) . 
Hence, we seek a proposal distribution that matches the CPHD 
cardinality exactly while exploiting the weights of individual 
labeled target clusters as computed from the approximate 
posterior PHD. A single component Vo-Vo density can be used 
for this purpose, 

q fe (X fc | X k -uz k ) = A(X fc )w(£(X fc )) [p(-)] Xt (52) 


Multi-Target Particle Filter 
with LMB Proposal Distribution 

Initialize particles ~ p 0 (-) 

For ft = 1 ,...,K 
• For i = 1,..., N p 

- For each t e £ (xjj.*E) 

* Generate a ~ ^[0,1] 

* If a < Pk,s(Z) g enerate ~ N(p-k,e, Qk,t ) 

- For each <e4 

* Generate a ~ ^[0,1] 

* If a < generate xfy ~ N {nk,i,Qk,t) 

- Evaluate the transition kernel 

(x«ix«i) = n (i-psW) n (i-psW)x 
n psw /4-i(4',ii4-i,() n 

(51) 


- Evaluate the proposal distribution 

qfcifc-^x^lxWj,^) = 

n n [}-pt B w)x 

Pk,s(^ X* (^ X k > Pk^iQk,^ x 

El Pfc.sW 

- Evaluate the multi-object likelihood g k 

- Update the particle weight w^' 1 using eq. jl9j 
• Normalize the weights and resample as usual 


B. Vo-Vo Proposal Distributions 

The LMB proposal distribution leads to an efficient sam¬ 
pling strategy for the multi-target particle filter. However, the 
LMB proposal does not exploit all the information from the 
SA-CPHD filter since the cardinality distribution of the LMB 
proposal does not match the cardinality distribution pk(ji) 
from the CPHD prediction/update. Matching of the cardinality 
distribution p k (n) is important if we are interested in designing 


We now specify the component weight w(£(Xfc)) and the 
multi-object exponential [p(-)] Xfe , needed to match the CPHD 
updated cardinality distribution and to account for the weights 
of individual target clusters. Clearly, the single-target densities 
p(-) are obtained straightforwardly from Gaussian clusters, i.e. 


Pk(x, t) = J\f(x; p k ,e, Qk,e), ( G L 0:fc (53) 


The weight w(£(Xfc)) is then chosen to preserve the CPHD 
cardinality distribution, and for a given cardinality, to sample 
labels proportionally to the product of the posterior PHD 
masses of any possible label combinations. Specifically, from 
the posterior PHD mass of each cluster p k (£) we construct 
approximate “existence” probabilities as 


r 


+ 

k 


U) 


PkU) 


3 = !>•••> |Lo:fc| 


(54) 


The cardinality of the set of labels Lo : fc, including birth 
and survival labels, grows exponentially in time. Moreover, 
in any practical implementation the use of a finite sample 
approximation coupled with resampling strategies typically 
leads to a much smaller unique labels set at each time k. 
Thus, eq. ( [54] ) is implemented by considering only labels from 
resampled particles at time k — 1, 


PkU) 


+ ( 7 ') 

' — |LoM +, 


EtTW) 


3 = 1 , • • • , |L 0 :lk 


Nr, 


L 0:fe =U^(X« 1 ) U L * 


(55) 

(56) 


The weight w(£(Xfc)) is then defined as 

r +f.)£(x fc ) 

u>(C(X k )) = Pk ( |£(X fc )|) fcU (57) 

e\c(TL k )\{ttk) 


where R k = {4"(j)}- e £ denotes the set of “existence” 
probabilities for all current tracks and e n (-) is the elementary 
symmetric function of order n. The construction of the pro¬ 
posal in ( |52| leads a simple and efficient strategy for sampling. 
Specifically, to sample from ([52]) we. 


sample the cardinality 


44 


of the newly proposed 


particle according to the distribution p k (ji ), 

sample xj^ labels £(xj^) from Lo : fc using the distri¬ 


bution defined by 


r+(-) £(Xk) 

e\c(x k )\(Rk)’ 












• for each t € £(Xj^) we sample the kinematic part x k \ 
from pk(-J) = N{",Vk,t,Qk,i)- 

A detailed pseudo-code for implementation is reported below. 

Multi-Target Particle Filter 
with Vo-Vo Proposal Distribution 

Initialize particles X*, 1 * ~ p 0 ( ■) 

For k — 1 ,..., K 

• For 1,..., N p 

- Sample the cardinality for the new particle |xj^ | ~ p k (n) 

- Sample the set of labels C j uniformly from L 0: fc 

- For each i e C (x^^ generate xfy ~ M (•; p, kj (, Q k ,e) 

- For j = 1_ ,N P evaluate the transition kernel 

ffcIfc—t (x«|X^r) = 

n a-psW) n (i-pbW)x 

lct, D K ' KzL ~' N B 

rr psw fk\k-i( x k ^\ x k ii : i ) x 

n pbw PB(?k\) 


that guarantees ffc|jfe_i(X^|X^ 1 ) ^ 0. This implies that the 
particle XjE has to verify the condition £(X^) nLo : fe_i C 
£(XjEi)> i.e. the labels set of the particle XjE includes the 

u\ 

labels set of surviving targets in the sampled particle Xjj, . 

V. Numerical Example 

In this section we demonstrate the RFS multi-target tracker 
with Vo-Vo proposal via a radar tracking application for 
closely spaced targets and low signal-to-noise ratio. 

A. Dynamic Model 

The kinematic part of the single-target labeled state vector 
Xfc = (xk,£k) at time k is described by x k = [xj),C s k\ T , 
which comprises the planar position and velocity vectors 
Xk = \px,k,Px,k,P v ,k,Py,k) T m 20 Cartesian coordinates, 
respectively, and the unknown modulus of the target complex 
amplitude g R. A Nearly Constant Velocity (NCV) model 
is used to describe the target dynamics, while a zero-mean 
Gaussian random walk is used to model the fluctuations in 
time of the target complex amplitude, i.e., 

x k +i = Fxk+Vk, v k ~ AT (0; Q) (60) 


- Evaluate the proposal distribution 

qfclfc-r (x^l X k _ u z k ) = 

w(£(X«)) nJV(xjj ) ;ttk,£,Qk,e) 

Pei 

- Evaluate the multi-object likelihood g k (z k |x£^ 

- Update the particle weight w jjy using 

N p 

g k (z k \x.% ) )^f klk - 1 (x.t ) \x.i i ± 1 )wl j ± 1 

w f =-- 

q fe (xW| X k _ u z k ) 

• Normalize the weights and resample as usual 


Notice from the update step in the pseudo-code that for 
each particle (z) we require the evaluation of the multi-target 
transition kernel with respect to the previous set of particles 


N„ 




fc- ! (Xfc I Xj._ 1 ) 


(58) 


3 = 1 


This is known as sum kernel problem ]47|-p9) and is due 
to the fact that the our proposal q/,. (Xp|X fe _!, Zk) does 
not depend on the previous particle XjE (as in the LMB 
proposal) but on the whole set of particles at k— 1. Efficient ap¬ 
proximation techniques exist for mitigating the computational 
load due to the sum kernel problem in ( |58j ), see g7j-gg. 
Furthermore, a simple approximate solution can be obtained by 
using a single particle xjE to evaluate the transition kernel 


X^-rflX^) ^-rflX^). (59) 

3=1 


where the matrices F and Q are defined as in |50|. 


B. Measurement likelihood function 

We now describe a multi-target observation model for multi¬ 
target tracking using radar measurements. A target x £ X 
illuminates a set of cells C(x), where C(x) is usually referred 
to as the target template. A radar positioned at the Cartesian 
origin collects a vector measurement z = [z^ ... z^ m l] 
consisting of the power signal returns, i.e. 




2 

E AWh$(x)+wW 




xGX:iGC(x) 


where z}^ is the complex signal in cell (z), and 

• is a zero-mean white circularly symmetric complex 
Gaussian noise with variance 2cr^ 

• /i^(x) is the point spread function in cell (z) from a 
target with state x 

fc«(x) = exp (_ (r‘-r(x)) a _ (* ~ d ( x ^ _ 

A w V 2 R 2D 2 B J 

( 62 ) 


where R, D and B are constants related to the radar cell 
resolution; r(x), d(x), and 6(x) are the target coordinates 
in the measurement space; and ri,di,bi are the cell 
centroids. 

• A(x) is the complex echo of target x, i.e. A(x) = 
A x e- ?s +a(x) with A x a known amplitude, 9 ~ [0 , 27t) an 
unknown phase, and a(x) a zero-mean complex Gaussian 
variable with variance a 1 , ,. 

a V x / 

For a non-fluctuating target amplitude (Swerling 0), A(x) is 
modeled as: 


However, in order to use the approximation in ( |59| we have 
to choose the index (to) of the previous particle in a way 


A(x) = A x e je , 9 ~ W [0) 2 tt) 


(63) 
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Let z ( ' v denote the deterministic part of the signal in cell i: 



Z {i) 

Z A 

2 

]T a x h { ;\ x ) 




xGX:zGC(x) 


The power measurement in cell (i) can be written as: 

2 b) _ |^W e is _|_ u ,| 2 

= (zj? cos (0) + 5R(w)^ + /za^ sin (6) + 7s(w)^ 

«= Ur + Uj 


where Ur ~ A/"(z^ cos(0), and Ui ~ sin(0), tr^) 

are statistically independent normal random variables. Then 
y/Ufi + Uj has a Ricean distribution, and reduces to a 
Rayleigh distribution when z/ = 0. Let SNR be the signal- 
to-noise ratio debited in dB as 


SNR = 10 



(64) 


We can choose crjj, = 1 so that A x = V2 io s,Ar ' R / 10 . In turn, 
since s/U^ + Uj ~ Rice(£^ , 1), the measurement in each 
cell ^(*) is described by a non-central chi-squared distribution 
with 2 degrees of freedom and non-centrality parameter 
and simplibes to an exponential distribution in the case = 
0. Then, the likelihood ratio for cell (i) is given by: 

<^(z^|X) = exp 0.52^ I 0 (^V (65) 

where lo(-) is the modibed Bessel function. Hence, the like¬ 
lihood function for the vector measurement z takes the form 


g(z|X) oc H <p(z® |X) (66) 

iec(x) 

where C(X) = C'(xi) U C(x 2 ) U ... U C(x|x|) is the union 
of all single-target templates, i.e. the set of measurement bins 
used for the measurement update. 


C. Simulation Results 

We consider a scenario with 3 incoming targets as depicted 
in Fig. [T] To better highlight the spacing of targets, we 
report in Fig. |T| the range-azimuth grid. Notice how the 
targets share cross-range cells for most of the simulation. 
Consequently, a correct estimation of the number of targets is 


very challenging. In Figs. 2(a) and 2(b) we report a snapshot 
of the linear domain power measurement in range-azimuth 
at time step k = 8 for both cases of SNR = 7dB and 
SNR = 10 dB, respectively. Relevant parameters used in 
simulation are reported in Table [I] For the biter initialization, 
we describe prior knowledge using a Gaussian distribution 
A I~( x b,Qb) which mainly proposes incoming particles (i.e., 
most of particle velocity vectors are directed towards the 
radar position). Specihcally, for the birth intensity we use 
the Gaussian Af (xg, Qb) and generate 5000 new born single¬ 
target particles in the CPHD biter at every time step. 

Results in terms of the estimated number of targets and 
Optimal Sub-Pattern Assignment (OSPA) distance m are 
reported in Figs. 3(a)|3(b) for the case with SNR = 10 dB 


Table I: Common Parameters 


Parameter 

Symbol 

Value 

Range Resolution 

R 

10 m 

Azimuth Resolution 

B 

i° 

Doppler Resolution 

D 

1 m/s 

Signal-to-Noise Ratio 

SNR 

(7,10 }dB 

Target Maximum Acceleration 

Umax 

1 m/s 

1 st Target Initial State 


[1260,-11,1240,-9] 

2 nd Target Initial State 


[1250,-10,1250,-10] 

3 rd Target Initial State 

x a 

x 0 

[1240,-9,1260,-11] 

Target Birth Mean 

X S 

[1250,-5,1250,-5] 

Target Birth Covariance 

Qb 

diag ([7.5,10, 7.5,10] 2 ^ 

Birth Probability 

Pb 

0.05 

Survival Probability 

Ps 

0.95 

n° of multi-target particles 

N v 

{3e3, 5e3,10e3, 20e3} 



x coordinate (m) 


Figure 1: Simulated scenario and estimated trajectories. The 
radar is positioned at the Cartesian origin and there are 3 
closely spaced targets moving towards the origin. 


and in Figs. 4(a)|4(b) for the more difficult case with SNR = 
7dB. Notice from Fig. 3(b) the increase in OSPA distance 
after k = 3 and k = 5, i.e. the instants at which new targets 
enter the scene, and then reduces in time thanks to the biter 
convergence. For the more difficult case with SNR = 7dB, 
we notice in Fig. 4(b) a slower biter convergence and in¬ 
creasing OSPA distance also for k = 15 and k = 20, 
i.e. the instants at which targets disappear. This is due to 
the fact that for lower SNR the tracker prefers to retain a 
“false” track for few steps rather then declaring a “dead” 
target too soon, as conhrmed by the time behaviour of the 
average estimated number of targets in Fig. |4(a) Tuning of 
the survival probability Ps can reduce this phenomenon. In 
practice, a perfect tuning of Pg and of the birth probability Pr 
require additional prior knowledge on the surveillance area. 
Overall, the results conbrm the applicability of the proposed 
approach for challenging multi-target problems with closely 
spaced targets and low SNR. 






































Estimated Number of Targets &T Estimated Number of Targets 


Azimuth Domain (deg) ™ Ra tt g e Domain (m) Azimuth Domain (deg)".= ' J 3p5i^'mT 1:| " jj^gj, D oma } n ( m ) 

(a) 7 dB for k = 8: Targets at (r = 1.7 km, b = 0.78 deg) (b) 10 dB for k = 8: Targets at (r = 1.7 km, b = 0.78 deg) 

Figure 2: Snapshot of the simulated radar power returns for low SNR and closely spaced targets. 
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(b) Monte Carlo results: OSPA Distance for SNR = 10d-B 


Figure 3: Monte Carlo results for the case SNR = 10 dB 




Figure 4: Monte Carlo results for the case SNR = 7 dB 



















































11 


VI. Conclusions and Future Research 

In this paper we discussed a general solution for multi¬ 
target tracking with superpositional measurements. The pro¬ 
posed approach aims at evaluating the multi-target Bayes 
filter using SMC methods. The critical enabling step was the 
definition of an efficient proposal distribution based on the 
Approximate CPHD filter for superpositional measurements. 
Numerical results confirmed the applicability to challenging 
multi-target tracking problems for closely spaced targets using 
radar measurements with low SNR. A large-scale application 
of this approach might not be possible due to worsening 
depletion problems in high dimensional state spaces. However, 
MCMC methods could be used to devise a particle implemen¬ 
tation that can scale with an increasing number of targets. 
Furthermore, subdividing the targets into statistically non¬ 
interacting clusters and then processing the clusters separately 
could lead to satisfactory performance with reduced computa¬ 
tional load. Finally, the capability of separating closely spaced 
targets for superpositional measurements, i.e. estimating the 
correct cardinality when there are unresolved targets, means 
the approach could also be used as an initialization block of 
cheaper trackers like the LMB |[32] and Vo-Vo f3l) biters. 
Specihcally, parts of the radar superpositional measurement 
could be processed with the proposed approach to bnd the 
correct number of targets as well as there location, then thresh- 
olded measurements could be processed with the LMB/Vo- 
Vo tracker. This should lead to improved performance as the 
LMB/Vo-Vo tracker would be using a more informative prior 
distribution. 
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